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Abstract. We study a connection between the isochronicity of a center of a 
polynomial vector field and the existence of a polynomial commuting system. 
We demonstrate an isochronous system of degree 4 which does not commute 
with any polynomial system. We prove that among the Newton polynomial 
systems only the Lienard and Abel systems may commute with transversal 
polynomial fields. We give a full and constructive description of centralizers of 
the Abel polynomial systems. We give new examples of isochronous systems. 



A center of a planar dynamical system is said to be isochronous if all cycles near it 
have the same period. The isochronicity problem has a long history After two or 
three decades of oblivion, the interest in this problem has been revived partially due to 
proliferation of powerful methods of computerized research. At this juncture, most attention 
is paid to polynomial vector fields (see, for instance, ^ and the bibliography therein). 
The characteristic feature of the present-day activity is the fact that some long forgotten 
results are rediscovered and at the same time new nontrivial examples of isochronous systems 
appear, giving rise to fresh and ingenious ideas. 

A Dynamical Meaning of the Cauchy— Riemann Conditions 

One of these ideas leans upon a remarkable observation by Villarini ||^ which allowed 
him to gave an elegant proof to a well-known theorem by Lukashevich [|l], ^ . Lukashevich's 
theorem asserts that each center of a system 

whose right-hand sides satisfy the Cauchy-Riemann conditions 

Px Qyy Py Qx 
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is isochronous. 

From the viewpoint of dynamics, these conditions mean that system (0) commutes with 
the orthogonal system 

y =-pix,y). ^ 

Indeed, the commutator, or Lie bracket of the vector fields (p, q) and (g, —p) is the field 
with the components 

QiPx-qy) -p{Py + qx), 
Pi.Px - Qy) + qiPy + Qx) 

whose vanishing at a noncritical point with + 7^ is equivalent to the fulfillment of the 
Cauchy-Riemann conditions. 

Thus, if a point, moving according to equations (|l|), returns to the initial position at time 
T and the orthogonal flow afterwards carries it to a nearby orbit of (|I|), then the fact 
that the flows commute implies that the motion of the point along the latter by time T also 
terminates at the same place where it starts. 

This argument proves Lukashevich's theorem and also allows us to make some interesting 
conclusions about the properties of the system under consideration. For example, the orbits 
of such a system, lying nearby a cycle, are necessarily closed, so that there are no limit 
cycles. Moreover, if an analytical system is defined over a simply connected domain then 
each cycle of the system surrounds a unique singular point, which is an isochronous center, 
and lies in the basin of this center together with the region that it bounds. 

We now consider an arbitrary system and suppose that it commutes with another 
system 

X = r(x, y), ~ 
y =s{x,y). ^'^ 

In the language of vector fields, this amounts to validity of the equalities 



rpx + spy - pr^ - qry = 
rqx + sqy - pSx - qSy = 



Then the above argument clearly remains valid in each domain where systems and (|^) are 
transversal. In particular, if such domain surrounds a center of system ([l|) then the center is 
isochronous. It was Sabatini ^ who noticed this fact. He also proved the converse assertion, 
namely: in a punctured neighborhood of an isochronous center there always exists a system 
transversal to the original system and commuting with it; moreover, this system has the 
same smoothness as the original system. This gives a new tool for proving isochronicity of 
various systems in a uniform way. 

It is worth noting that conditions amount merely to local commutation of the flows 
Ut and corresponding to systems (|l]) and @, even if the corresponding vector fields are 
defined on the whole plane. Surely, sometimes we can extend equality between the values 
VaUr and UrVa at a given point to rather large r and a, for example, when the products 
V^'Ur' and Ur'V(j> make sense at this point for all pairs o"', r' in the rectangle < a' < a, 
< t' < T. However, we cannot guarantee the equality Va-Ur = UrVa- on the only reason that 
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both sides of the equahty are defined at the point in question. Neglect of this circumstance 
might lead to erroneous conclusions. For example, by an argument similar to that above, 
one might conclude that two isochronous centers whose basins are joined by an orbit of a 
commuting system have the same period. However, this is not always the case. We exhibit 
one example. 

A system z = f{z) on the plane, with f{z) a holomorphic function of a complex variable 
z, has a center at a point zq if and only if /(-Zo) = and the derivative f'{zo) is a nonzero 
purely imaginary number. This result is also due to Lukashevich f^, but essentially the 
same result was recently reproved in 0. Thus, if, say, f{z) = iz{l — z^) then our system 
has three centers. They all lie on the real axis, separating it into four parts each of which, 
as we easily see, serves as an orbit of the orthogonal system z = —if{z). Two of them join 
the point (0,0) to the points (—1,0) and (1,0). Moreover, the extreme centers have period 
TT, whereas the period of the middle center equals 27r. 

We now study the mechanism of this phenomenon in more detail. The essential elements 
of the phase portrait of the system under consideration are displayed in Fig. 1. Two sep- 
aratrices infinite to both sides separate the plane into three domains. These domains are 
entirely covered by cycles. The whole picture is symmetric about the coordinate axes. 




Fig. J Fig. 2 

In the orthogonal system, consider the motion starting at t = at the point A indicated 
in the picture and finishing at some time t = a at the point B lying at the intersection of one 
of the orbits around the right center with the abscissa axis. If we use the former notations 
U and V for the original and orthogonal fiows and take symmetry of the system and the 
above- indicated values of periods into account, then we come to the following expressive 
relations: 

UMA) = U^{B) = B, V^UM) = V.{-A) = -B. 

Certainly, the culprits are the separatrices separating the basins of the centers: their 
blame consists in their going to infinity too fast, so that they create obstacles to extending 
commutation of the fiows to the needed range. Indeed, we can easily calculate that these 
separatrices represent two branches of the hyperbola 2x'^ — 2y'^ = 1 and the motion along 
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them is described by the exphcit formulas 



cos(t/2) sin(t/2) 
X = ± . =, y = ± . 

V 2 cost V 2 cost 

The duration of traversing each of these orbits embeds in the interval — 7r/2 < t < +7r/2 
of length IT. Thus, if a* is the moment when the point C = Vcr^{A) occurs on the right 
separatrix then the orbit Ut{C) terminates by time t = 7t/2. Now, disregarding the fact 
that the other orbits of the flow U which start at the points of the segment AB live for the 
infinite time, there is no "commuting" rectangle which could equate the values of VcrU-,,- and 
U^V^ at A (Fig. 2). 

Seizing an opportunity, we indicate one more curious property of the periods of polyno- 
mial analytical systems. 

Theorem 1 If a polynomial system of degree n >2 satisfies the Cauchy-Riemann conditions 
and has n singular points such that n — 1 of them are centers, then the remaining point must 
he a center too. Furthermore, among all singular points there necessarily exist those around 
which the rotation is clockwise as well as those around which the rotation is counterclockwise. 
Finally, the sums of periods for the various types of centers are the same. 

Proof. If among the roots zi, . . . , z„ of a polynomial f{z) of degree n > 2 there are 
n — 1 simple roots then it is clear that the remaining root is simple as well. Furthermore, 
the following remarkable identity is easy to verify: 



If all points Z}. with k = 1, ... ,n — 1 are centers of the system z = f{z) then, as was 
mentioned above, the numbers f'{zk) are purely imaginary. But then the preceding identity 
implies that the derivative f'{zn) must be purely imaginary as well. Since this derivative 
differs from zero, the point z„ must be a center. It remains to observe that the absolute value 
of the real number 2TTi/f'{zk) equals the period of cycles around the (isochronous) center Zk 
and the sign of this number determines the direction of the circuit. 

Sabatini's Problem and Devlin's Example 

We now discuss the application of the above approach to studying isochronous centers of 
polynomial dynamical systems. The following problem was stated in the already-quoted 
article 0: 

Problem. Is it true that a polynomial system, having an isochronous center, commutes 
with some transversal polynomial system! 

For every quadratic system with an isochronous center, in the author of this question 
constructed a commuting polynomial system whose linear part is orthogonal to the original 
linear system. Using Mathematica [|1^, we examined all isochronous polynomial systems 
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described in the articles we know. It is surprising, but for the majority of them the answer 
to the above question turned out to be affirmative. We outhne the results of our inspection. 

Example 1. The article contains a full list of all (up to homothety) cubic systems 
symmetric with respect to an isochronous center. It is as follows: 

X = —y — Sx'^y + y^, x = —y + x'^y, 

y =x + x^ — 3xy'^, y =x + xy'^, 

X = —y + Sx'^y, X = —y — Sx'^y, 

y = X — 2x^ + dxy"^, y = x + 2x^ — 9xy'^. 

Below we indicate the corresponding polynomial systems, obviously transversal to the 
original systems at least near zero, which commute with them as one can check by straight- 
forward verification of identities (^): 

/y» <-y» I <-v»3 '^'T''?/^ If If /y»^ 

y =y + Sx^y -y^, y =y- x'^y, 

X = X — 5x^ + 6x^, X = X + 5x^ + 6x^, 

y = y — 9x^y + ISx^y, y = y + 9x^y + ISx'^y. 

Example 2. As is well known among Kukles's systems we have exactly one system 
(up to a similarity transformation) with an isochronous center: 

X = -y, 

y = X + 3xy + x^. 
For it, there exists a transversal commuting polynomial system: 

X = X + xy + x^, 

2 I 2 4 

y = y — X + y — X . 

As computer analysis shows, the above isochronous system remains unique even in the full 
class of Kukles's systems when the second equation is specified by an arbitrary polynomial 
of third degree. Some relevant remarks can be found, for instance, in 

It is curious to notice that each of the isochronous systems of the shape 

X = -y +p^[x,y), 
y = x + q2ix,y), 

where p2 and q2 are homogeneous quadratic polynomials written in the form they are pre- 
sented in Table 1 of 0, in the variables X = x, Y = y —p2{x, y) coincides, up to a similarly 
transformation, with the above-written Kukles's system. 

Truly ubiquitous, this system plays a magical role in the problem under study. We 
encounter it once more in our article, moreover, at a place where we did not expect to see it 
at all. 
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Example 3. Studying the cubic Kolmogorov systems, the authors of [|TT1 arrived at the 
following system with an isochronous center at the origin: 



X = —y + 2xy — ax'^y, 
y = X — x^ + y"^ — axy^. 



This system commutes, for example, with the polynomial system 

X = X — x^ + y"^ — axy^, 



y = y — 2xy — ay^. 

It is noteworthy that, in all three examples of systems, their direction fields are sym- 
metric about the abscissa axis, which makes it obvious that the origin is a center of each 
of them. Having constructed transversal commuting systems, we obtain one more proof for 
isochronicity of these centers. 

Example 4. A new family of isochronous cubic systems was found in 0. Slightly 
changing notations, we can write it as follows: 

X = —y — 2xy + x^{a + ax + Py), 
y = X + x^ — y"^ + xy{a + ax + (3y). 

For arbitrary real a and /3, the origin is a center of the system, which is easy to check 
by using the first integrals indicated in 0]. Incidentally, everything happened just the other 
way about: the authors started with choosing an integral and then wrote a system for it. 
They studied the result by the Darboux linearization method and proved isochronicity of 
the center. Another proof of isochronicity appears on observing that the system in question 
commutes with the following: 

X = X + x"^ — y"^ + xy{a + ax + Py), 
y = y + 2xy + y'^{a + ax + [3y). 

As we can see, this example contains the previous: it suffices to take a = and (3 = —a 
and change the sign of the variables x and y. 

Example 5. A remarkable class of isochronous hamiltonian systems relates to area- 
preserving transformations |T^. Suppose that functions u = u{x,y) and v = v{x,y) satisfy 



the identity 



UxVy - UyVx = 1. 



In that event, if m = f = 0, say, at x = y = then the origin is an isochronous center of the 
system 



X = —UUy — VVy, 

y = uu^ + vv^. 

This is obvious, since in the coordinates u, v the system takes a standard form of the 
harmonic oscillator: u = —v, v = u. However, for us it is important to indicate an explicit 
form of a commuting system which is polynomial in case u and v are polynomials: 

X = UVy — VUy, 

y = -uv^ + vu^. 
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The fact, that this system really commutes with and is transversal to the original system, can 
be verified by straightforward calculations, but it is simpler to observe that in the coordinates 
u, V the system is written as the system ii = u, i) = v corresponding to the canonical Euler 
field {u, v) obviously orthogonal to the field {—v, u) of the harmonic oscillator and commuting 
with it. 

It is easy to construct as many isochronous hamiltonian systems as we wish if we observe 
that the Jacobian of the mapping (m, f ) is preserved under the replacement oiu by u + ip{v) 
or w by f + V'('w); where ip and if) are arbitrary smooth functions of one variable. 

We thus see that the hope of a positive solution to Sabatini's problem had certain grounds. 
However, the existence of a commuting system, even in a small neighborhood of a singular 
point, imposes strong requirements on the global behavior of the system as soon as we talk 
about the polynomial vector fields. It is quite possible that these requirements are met by 
all cubic isochronous systems, but they are already incompatible with the degree of freedom 
of systems of forth degree. One of the circumstances in support of what was said is refiected 
in the following simple theorem slightly developing the arguments of the beginning of the 
article. 



Theorem 2 // two polynomial systems commute in some domain and are not collinear in 
at least one point of the plane then every center of each of the systems is isochronous. 

Proof. The points where a polynomial of two variables vanishes either cover the whole 
plane or constitute a closed nowhere dense set. Therefore, first, our systems commute on all 
of the plane and, second, the basin of every center of each of them includes a cycle that the 
other fiow, commuting with the given system, intersects somewhere transversally. But then 
all familiar arguments apply and demonstrate that all nearby cycles have the same period. 
In this case the period, depending analytically on the orbit, has a constant value on the 
whole basin. 

Thus, the natural question arises: Can two centers of a polynomial system coexist if one 
of them is isochronous and the other is not? In view of the preceding theorem, a positive 
answer to this question would be simultaneously a negative answer to Sabatini's problem. 
We found an appropriate example in the article WM, wherein the system 



X = —y — x"^ + 4x^2/^ + y^, 
y = X — 4x^y 

is indicated which has one isochronous and two nonisochronous centers. 

When our article was in print, we received from Sabatini a large, still unpublished survey 
of isochronous centers [0] which contains many of the commuting polynomial systems we 



consider, as well as a whole series of new examples. Interestingly, this survey also points to 
Devlin's example []T3|, moreover, in connection with the same question. 
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Homogeneous Perturbations of the Harmonic Oscillator 



Besides the above, there are other reasons for which an isochronous polynomial system may 
fail to commute with any nonproportional polynomial system. Each of the following systems 
can serve as an example: 

where m is an arbitrary natural number. Here we have only one singular point and therefore 
the preceding arguments fail. 

We have found and investigated this system before we learned about Sabatini's problem. 
Now, when the question has received an answer, we might confine ourselves merely to men- 
tioning our example. However, since polynomial vector fields of degree higher than tree are 
not studied practically and any information could be useful, in this section we conduct a 
small investigation of system which is curious in many respects. 

Theorem 3 The origin is an isochronous center of system . 

Proof. Symmetry of the direction field of {ref5) about the ordinate axis implies that 
the origin is a center of the system. We write the system in polar coordinates: 



g = Q^^+' cos'""-' ^ sin' ^, ^ = 1 + 0'""+' cos'"" ^ sin ^. 
Differentiating these equations with respect to t, we can easily verify that the function 



t-^- 



2m + 1 



remains constant along each orbit. If the orbit is closed and, in consequence, surrounds 
the origin, then during one circuit the angle "i? receives the increment 27r and the fractional 
summand of the above expression returns to the initial value. Therefore, the duration of the 
circuit equals 27i. 

Theorem 4 The polynomial systems able to commute with system @ may differ from 
only by a constant factor. 

Proof. Suppose that the field {p, q) of system (5) commutes with a polynomial vector 
field (r, s). When p = —y, as we have in our case, the first expression in (4) takes the shape 

s = yr^- qry. (D 

Now, if we rewrite the second relation in @) by using (|]) and next compare summands from 
the viewpoint of their degree in y, not forgetting that the degree of the polynomial q is 
exactly three, then we find that the degree r cannot exceed one, so that 

r = ro(a;) + ri{x)y. 
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The further, rather elementary examination leads to the equalities tq = and ri = — c, 
where c is a constant. Hence, r = —cy and so s = cq. Thus, the field (r, s) is coUinear to 
the field {p, q) with a constant coefficient of proportionality. 




Fig. 3 

In Fig. 3 we depict the scheme of the phase portrait of system (j^) which is obtained by 
the classical methods of the qualitative theory of planar dynamical systems. The above ex- 
ample of an isochronous system resulted from studying the one- dimensional nonconservative 
Newton equations with homogeneous nonlinearities added to the equation of the standard 
harmonic oscillator. On the phase plane, these equations are expressed by the system 

y =x + Qn{x,y), ® 

where Qn{x,y) is a homogeneous polynomial of degree n. Examining the known full lists 
|]T|, ^ of isochronous systems involving only quadratic or only cubic nonlinearities, we can 
observe that among them there are no Newton systems of the form with n = 2 and n = 3. 
Inspecting the necessary conditions for isochronicity of a center by means of Mathematica, 
we verified that, up to homothety, there is exactly one isochronous system for n = 4: it is 
system (^) with the exponent m = 1, and if n = 5 then the origin never can be an isochronous 
center of system (^. By this we once more confirmed the results that were obtained thirty 
years ago in 
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In the more difficult case of n = 6, we succeeded only in examining systems that are 
symmetric about one of the coordinate axes. We found out that, with the previous stipu- 
lations, the only system among them which has an isochronous center at zero is system (^) 
with m = 2. Finally, if n = 7 then system cannot have an isochronous center at zero. 

The above consideration and some other reasons allow us to pose the following two 
questions. 

Problem 1. Is it true that a system of the form with nonlinearities of degree higher 
than two has a center at the origin if and only if its direction field is symmetric about one 
of the coordinate axes'? 

Problem 2. Is it true that the origin is an isochronous center of system (0) only if the 
system has an even degree and is reduced to system by a similarity transformation? 

Surely, in near future the progress in computers will advance the study of the class of 
systems in question far beyond the point we have reached. However, in this way we could 
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hope to receive only a negative answer to the posed questions. If the answers are affirmative 
then the attempts at their justification will probably lead to new profound ideas. These 
questions are also closely connected with the other problems to be dealt with below. 

Polynomial Newton Systems 

The circumstances mentioned in the proof of the last theorem gives us the thought of try- 
ing to describe all polynomial Newton systems for which there are transversal commuting 
polynomial fields. The results turned out to be unexpected for us. The rest of the article is 
devoted to exposition of these results. 

Consider a general polynomial Newton system, expanding its force function in the powers 
of the variable having the meaning of velocity: 

^ " (|8D 

ij = qQ{x) + qi{x)y + . . . + qn{x)y'^ . ^ 

The polynomials . . . I'n. x can have arbitrary degrees, but we assume that go and qi 
have no constant terms and that the derivative go zero equals one. This means that the 
linear part of the system is the canonical oscillator. 

As is well known, if a polynomial Newton system (|^) is conservative, i.e., all g^ with k > 1 
vanish; then the oscillations described by the system are isochronous only in the trivial case 
of qo{x) = X. This fact was ffist established in |TB|, afterwards it was indicated in ||T|, and 



after many years it appeared in [|T7|, where it received a new proof. Thus, in the conservative 
case, there arises no question of existence of a commuting system. Therefore, henceforth we 
assume that n > 1 and that the last polynomial g„ does not vanish identically. 

Theorem 5 If n > 4 then system (§) cannot commute with any polynomial system nonpro- 
portional to it. 

Proof. Suppose that a polynomial vector field (r, s) commutes with the field {—y,q), 
where g is the right-hand side of the second equation in (j^). As mentioned, the first relation 
in (^) allows us to express s by means of (^. In this case the second relation in (|D takes 
the form 



qx{r - yry) + {yqy - q)r^ + y'^r^^ - 2yqr^y + q^Vyy = 0. ® 

As in the proof of Theorem 4, this equality implies that the degree of the polynomial r 
in the variable y cannot exceed one, so that r = ro(x) + ri{x)y. Thus, the last summand 
in d^) vanishes and leading is the coefficient of ?/"+^ equal to (n — 3)g„r^. Since n > 3 and 
the polynomial g„ is not equal to zero, we conclude that r[ = or ri = — c, where c is a 
constant. 

Now, the maximal power of y in (^) equals n and the corresponding coefficient equals 

g^ro + {n- l)qnrQ. 
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Its vanishing is possible only if (ijn'^o"^)' — Oj hence, tq too is a constant. It remains to 
observe that in these circumstances equation (D is reduced to the form q^VQ = 0, where 
= I for X = y = 0, which implies that tq = 0. Thus, r = —cy, and so s = cq, as it follows 
from (^. 

Problem 3. Can a polynomial system of the form (H) with n > A be isochronous^ 

Anyway, such systems are absent from the list of the examples described in the articles 
known to us. 

Theorem 6 If n = 2 then a system of the form (H) commutes with no polynomial system 
nonproportional to it. 

Proof. Here we use the same notations as in the proof of the preceding theorem. Observe 
that the highest power with which the variable y enters in each summand of for n = 2 
equals m + 2, where m is the degree of r in y. Hence, 

(1 - m)q'2rm + (1 - 2m)g2?^m + C + ^{"^ - ^)(ll'^m = 0, 

with rjn the coefficient of in the expansion of r, representing a polynomial in x. Since q2 
is not identically zero, the last equality implies that m is at most one. Certainly, it suffices 
to settle the case of m = 1. Putting as before r = ro(x) + ri{x)y, we note that the equation 
is then transformed into r" = q2r[, which can hold only for a constant polynomial ri. 

Collecting the coefficients of in (^, we now arrive at the relation Tq = —{q2ro)' which 
means that ro too is a constant. The further arguments are reduced to repeating what was 
said in the proof of the previous theorem: rg = and the vector field (r, s) is proportional 
to the field {—y, q) of the system under study. 

Problem 4. Are there isochronous polynomial systems of the form for n = 21 

We hope that answers to the questions posed in this section will be related to a rather 
deep penetration into the isochronicity nature of polynomial systems. But already now it 
should be emphasized that the very statement of these questions became possible due to the 
remarkable idea that had combined concepts that, as it would seem, have a little in common: 
isochronicity of oscillations and commutation of vector fields. 

Isochronous Lienard Systems 

Of the two remaining classes of polynomial Newton systems, we first discuss one for which 
n = 1. In this case the force function depends linearly on velocity and the system, known 
under the name of a Lienard system, acquires an especially attractive simple form: 

X = -y, 

y = qoix) + qi{x)y 

We recall that the linear parts of the systems under study represent the standard harmonic 
oscillator. 
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Here we digress slightly from the main theme of our investigation, since in our opinion 
the question of commuting fields for systems of the form ([1^) is not very interesting for the 



following reasons: In Sabatini studied Lienard systems with odd analytical coefficients Qq 
and qi and obtained an exhaustive solution of the isochronicity problem in this case. Namely, 
he established that under the indicated conditions the origin is an isochronous center if and 
only if 

P(x) } 

qoix) = X -\ where lix) = / sqi(s)ds. 

x-^ J 



Furthermore, there are rather solid grounds to believe that, as far as the polynomial Lienard 
systems are concerned, there are no other cases of isochronicity, although nobody can prove 



this yet (in this connection, see the survey \T^] and the references therein). Anyway, we 



are firmly convinced that, among all polynomial Lienard systems, only those belonging to 
the Sabatini class can commute with some transversal polynomial systems; furthermore, 
they comprise a small part even in this class. Clearly, in such situation, the question of 
commutation looses actuality, unless its study finds the end in itself. 

We now formulate two problems whose positive solution would clearly speak in favor 
of the above-mentioned conjecture. Assume as before that go and qi are polynomials. We 
arrange their primitives: 



X X 

Qo{x) = J qo{s) ds, Qi{x) = J qi{s) ds. 



Problem 5. Is it true that the origin is a center of system (|r^) only in the following two 



? 



cases: either Qq and Qi are even polynomials or Qi is a polynomial of Qq 

This problem is cast by a beautiful result by Cherkas that he obtained in [[1^]. His 
theorem, on the one hand, guarantees the presence of a center under each of the above 
conditions and, on the other hand, enables us to reduce the problem to the following question: 
is it true that if a polynomial of one variable is represented as the sum of a power series in 
another polynomial which starts with a square of the variable and involves at least one odd 
power of this variable then the series terminates somewhere? 

Problem 6. Is it true that if Qi is a polynomial of Qo then system cannot be 
isochronous^ 

As an illustration of the second question, we consider the simplest class of systems sat- 
isfying the indicated condition: 

X = — I/, 

y =f{x){l + y), 

where f{x) is a polynomial starting with x. We prove that the origin is always a center of 
the system, but it is never isochronous. 
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X 



A' 



Fig. 4 



Fig. 5 



First of all, by setting 



X 



U{x) = J f{s)ds 



T{y) =y - \n{l + y), 







we find that the sum 



U{x)+T{y) 



-{x^ + y'') + ... 



serves as a first integral of our system. Next, we should use the fact that the function T(?/), 
defined for y > —1, has two intervals of strict monotonicity (Fig. 4) on each of which it takes 
all values from to +oo, and examine two possible cases: whether or not the system has 
stationary points other than zero. 

In the first case, until > attains the value of the polynomial U at the critical point 
of U nearest to zero (Fig. 5, a), the equation U{x) + T{y) = E determines two functions y 
of X whose graphs form together a cycle of the system. This means that the singular point 
nearest to the origin lies on the boundary of the basin of the center and so this center cannot 
be isochronous. If there are no stationary points other than zero, then for the same reasons 
the basin of the center (Fig. 5, b) adheres immediately to the line y = —1 representing 
an orbit of the system along which the motion has unit velocity and, consequently, infinite 
duration, so that a long time is needed for traversing the nearby cycles. 

The assumption that f{x) is a polynomial is needed only for the primitive U{x) to have 
two branches going to infinity in the case of the absence of nonzero roots and also for the 
period of a cycle around a center to depend analytically on the cycle. 



We now consider the polynomial Abel systems, more precisely systems of the form (§) with 
n = 3, slightly changing their record: 

X = -y, /p-^ 

y = qo{x) + 3qi{x)y + q2{x)y'^ + q3{x)y^ . 

We do not know all isochronous Abel systems, but we can give a full description for those 
of them for which there are transversal commuting polynomial fields. 



Centralizers of the Abel Systems 
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Theorem 7 The field of system (0) commutes with some transversal polynomial vector 
field if and only if the following two conditions are satisfied: 



qiqo = 3gf -q'o + l, 

q^ql = qf ~ q'lQo + qi- 



Proof. Let (r, s) be a polynomial field commuting with the field (— |/, q) of system 
([TT|). Then s, as was constantly mentioned, is expressed by formula (P) and r is found from 
equation (^. Repeating word for word the arguments presented at the analogous place in 
the proof of Theorem 4, incidentally, also devoted to one of the equations in the Abel class, 
we conclude that in this case as well we have r = Tq + ri?/, where Tq and ri are polynomials 
of X. 

Elementary calculations of the coefficients of various powers of y in (|^) lead us to four 
relations, the first two of which are as follows: 

q'^TQ - qoTQ = 0, q[ro - qor[ = 0. 

The general solution of these equation have the form 

'^o = cogo, n = cogi + ci, 

where cq and ci are arbitrary constants. If cq = then it is easy to see that the commuting 
field is proportional to the original field. Therefore, henceforth we may assume that cq = 1 
and ci = 0. 

Now, if we write down the third of the above-mentioned four relations, 

q2ro + 92^0 + - 6qir[ = 0, 
by using the equalities tq = go and ri = gi, we obtain the equation 

gko + q2qo + q'o - Qqiq[ = (g2go)' + (gd)' - ^{qiY = o 

which is easy to solve and which leads to the formula 

g2go + go - Ml = C2, 

where C2 is a constant. Factually, C2 = 1, since, as we have agreed, go(0) = gi(0) = and 
q'oi^) — 1- This proves the necessity of the first condition in (|T^. 

To prove the second condition, we use the last relation which ensues from (|]): 

q'^ro + 2g3ro + r" - q2r[ = 0. 

Expressing ro, ri and g2 in terms of go and gi, we obtain the differential equation 

qko + 2g3gogo + gog" + gogi - 3qfq[ -q[ = o 
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which integrates without effort: 

ml + -qf-qi = c3. 

Here C3 is one more constant. But now it relates to an expression vanishing at x = 0, so 
that C3 = 0, and the necessity of the second condition in ( ]T^ ) is estabhshed as well. 

Sufficiency of (|1^) for existence of a sought commuting system follows from the obvious 
conversability of our logical constructions. Also, it can be easily verified by inserting the 
found r in (|^). 

As we see from the proof of the theorem, in the case under consideration there is exactly 
one sought commuting field (r, s) which is, up to a constant factor and a trivial summand, 
proportional to the original Abel field {—y, q)'- 

r = qo + qiy, s = yr^ - qvy. (|TB|) 

Straightforward application of Theorem 7 to construction of commuting polynomial sys- 
tems encounters certain difficulties: if we choose the polynomials go and qi at random then 
the functions q2 and gs calculated in accordance with formulas ([T^ ) turn out as a rule merely 
rational. Therefore, it is desirable to have an explicit description for all polynomial Abel 
systems whose coefficients are bound by relations (0). This problem turned out to be 
solvable. 



Theorem 8 A polynomial Abel system (|TT]) commutes with some transversal polynomial 
system if and only if it has the following form: 

^ ^ (|TD 

y = {x + a'^x^){l + h{x)y)^ + 3axy{l + h{x)yy — h'{x)y^, ^-^ 

where a is an arbitrary number and h{x) is an arbitrary polynomial. 

It is easy to see that the coefficient of y^ in (|T^ vanishes only if h{x) = in which case 
the system transforms into the isochronous Kukles system. Even in the simplest case in 
which a = and h{x) = 1, we obtain a new example of an isochronous Abel system: 

X = -y, 

y =x{l + y)^. 

In this case it is obvious that the origin is a center. Isochronicity of the center is guaranteed 
by the membership of the system in the class (0), although here this fact can be verified by 
straightforward calculations. 

We turn to proving Theorem 8. In all forthcoming lemmas, the talk is about the polyno- 
mial coefficients of the second equation of the Abel system (pH]) which satisfy relations ([T2|) 
of Theorem 7. 



Lemma 1 The polynomial go has the form qo{x) = x + Ax^, where A is a constant. 
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Proof. The second condition in ([121) shows that the polynomial gi(l + qf) is divisible 
by go- Then it is easy to see that go admits the decomposition go = piP2 into the product of 
the polynomials 

= X + O(x^) and p2(x) = 1 + 0(x) 
which divide gi and 1 + gi respectively. Since 

g2go = 3gi - go + 1 = 3{ql + 1) - go - 2, 

the difference gp — 1 is divisible by pi and the sum + 2, by p2- Using the equality 
Qo = P1P2 + Pip'21 we arrive at the following important conclusion: there are polynomials hi 
and h2 such that 

p\p2 - 1 = hipi, P1P2 + 2 = h2P2- dUD 

Before going on further, we exclude one particular case. Namely, as the first equality in 
([TSp shows, if P2 = 1 then p'^ — 1 = hipi, which is possible only for pi = x; whence, go = x. 
This is in a full agreement of the claim of the lemma, if we take A = 0. Thus, the degree 
of p2 can be assumed positive. But then, as follows from (0), the degrees of hi and /12 are 
less exactly by one than those of p2 and pi. 

We now multiply the first relation in (p!5| ) by h2 and express the resultant product h2P2 by 
using the second relation. Interchanging the words "one" and "two" in the last proposition, 
we perform an analogous procedure again. At a result, we obtain two interesting equalities: 

2p[ -h2= Pi{hih2 - p'iP2), -p'2 + 2hi = P2{hih2 - p'lP^)- 

They are interesting, because they make their left-hand sides divisible by polynomials of 
higher degrees and so yield the formulas 

2hi = p^, h2 = 2p[. 

In this case, both relations of (|T5|) become a sole equality: 

2p[p2-pip'2 = 2. 

It is from the last relation that we derive the necessary conclusions. 

First of all, for the mutual compensation of the coefficients of the highest powers of the 



variable on the left-hand side of (16), it is necessary that the degree of p2, being nonzero, 
must be twice as large than the degree of pi. In that event, we can choose a factor A so that 
the degree of Ps = P2 — Ap\ be strictly less than that of p2- However, equality (|1^) remains 
valid if we substitute p^ for p2 in it: 

ViPa - = 2. 

These arguments about degrees now lead us to another conclusion: p^ = c, where c is a 
constant. Thus, p2 = Apf + c, and (16) is extremely simplified: cp[ = 1. Involving the initial 
data, whence we infer that c = 1. Therefore, pi = x and p2 = Ax'^ + 1, which completes the 
proof of the lemma. 
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Lemma 2 There must exist a real number a and a polynomial h{x) such that A = a? , where 
A is the coefficient mentioned in the preceding lemma and qi has the form 

qi{x) = ax + {x + a^x^)h(x). 

Proof. In view of Lemma 1, the first relation in (|12|) looks like 

q2{x + Ax^) = 3{ql - Ax^). 

Were the coefficient A negative, the left-hand side of the equality would have nonzero roots, 
whereas the right-hand side would not. This proves the possibility of writing A = a"^. Then 
the preceding relation can be written otherwise: 

q2x{l + a^x^) = 3(gi — ax){qi + ax). 

It is clear that one of the nonconstant factors on the right-hand side must be divisible by the 
irreducible polynomial 1 + Since the sign of the number a was not fixed, we may assume 
that this is the first factor. Moreover, the chosen factor is also divisible by x. Thus, the 
function h{x), necessarily defined as the quotient of qi{x) — ax hj x + a'^x^, is a polynomial. 

The next lemma, whose proof is reduced to simple calculations, is stated only for the 
reader's convenience. 

Lemma 3 // go (^i^'d qi are expressed in terms of a and h as it is said in the previous two 
lemmas, then equations (12) is uniquely solvable for q2 and q^ in the class of polynomials 
and their solutions have the form 

q2{x) = 6axh{x) + {x + a'^x^)h'^{x), 

qslx) = 3axh'^{x) + {x + a'^x^)h^{x) — h'{x). 

To complete the proof of Theorem 8, it now remains only a little: we have to insert the 
expressions for coefficients which are indicated in the three lemmas into the second equation 
of (|ll|) and verify that the result is system (|^) . Also we should observe that these coefficients 
satisfy relations (|T^), because it is from exactly these relations that the last two of them 
were found. 

Conditions for a Center in Smooth Abel Systems 

Describing polynomial Abel systems commuting with transversal polynomial fields, we sup- 
posed nothing about their behavior near the origin in advance. Meanwhile, if the linear 
part of a system is a harmonic oscillator then the origin can be not only a center but also a 
focus. In the latter case, it is quite possible that there exists a nontrivial commuting field, 
in which case we also observe some kind of isochronicity, namely: a point moving along a 
spiral around a focus returns to an arbitrarily given radial orbit of the commuting system 
by the same time. 

Thus, we have to clarify which of the found Abel systems have a center at the origin. 
Furthermore, it is interesting to understand the meaning of the enigmatic form of these 
systems which is indicated in Theorem 8. Answers to these questions are given in the 
following theorem which holds not only for polynomial systems. 
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Theorem 9 For an arbitrary number a and an arbitrary smooth function h{x), the system 
of the form (I4) in the variables 

X = x. Y 



1 + h{x)y 

is collinear to the isochronous Kukles system: 

XX = -Y, 

\Y = X + 3aXY + a^X^ 



where A = 1 — h{X)Y . In particular, if a = then the system is simply reduced to the 
harmonic oscillator. Moreover, the origin is always an isochronous center of system . 



Proof. Verification of reducibility of one system to another when a change of variables 
is given exphcitly requires nothing but accurate calculations. CoUinearity of a system to the 
Kukles system, although with a variable factor, shows that the origin is a center. Isochronic- 



ity of the center follows from the fact that the field of system (14) commutes with the 
transversal field (r, s) if we define the latter by formulas (p!3|). This is easy to check by direct 
calculation of the Poisson bracket or, equivalently, by verification of equation (P). 

The just proven theorem can be seen as an assertion refiecting a curious property of 
the isochronous Kukles system. Indeed, if we multiply or divide the field of this system 
by an arbitrary function not taking zero values in a neighborhood of the center then the 
orbits of the new system are certainly the same but the law of motion along them changes. 
However, as we see, if this field is divided by a function of the form 1 — h{X)Y, where h{X) 
can be arbitrary, then the total time of circuit around the center remains the same and 
isochronicity is preserved. In our opinion, this extraordinary phenomenon which, judging by 
all appearances, should also be encountered in other systems deserves an independent study. 
To state it more precisely, it would be interesting to clarify the changes of the velocity of the 
motion along orbits of an isochronous system which do not infiuence the period of cycles. 

Formulas (|T^) reduce system (|ll|) to the form (|14D only in the case of polynomial vector 
fields. If we remove this restriction and consider the Abel systems with arbitrary smooth 
coefficients then conditions (0) distinguish an essentially wider class of systems. It turns out 
that in the general case as well these conditions not only guarantee existence of a transversal 
commuting field but also ensure the presence of a center at the origin. 

Theorem 10 For arbitrary smooth functions qQ{x) = a; + 0(a;^) and qi{x) = 0{x), relations 
([T^ ) uniquely determine smooth functions q2{x) and q^lx) which are also analytical if so are 
qo{x) and qi{x). Moreover, each Abel system (11) with these coefficients has an isochronous 
center at the origin. 



Proof. Since the order of qo{x) in x is exactly one; to prove solvability of equations ([T2D 
in the class of smooth or analytical functions, it suffices to observe that the orders of their 
right-hand sides are at least 1 and 2 respectively. 



18 



It is curious to notice that systems (0) with coefficients connected by relations (|12|) are 
contained in the class of integrable Abel systems which was recently found by Lukashevich 



and described by him in [^. Following his constructions, we can obtain a ffist integral of 
our system which turns out to be equal to the sum of the squares of two functions X and 
Y of the variables x and y. Omitting the details of intermediate calculations, we only state 
the result. Namely, we define these functions by the formulas 

X = xe^W, Y = ^^y^ , 

qo{x) + qi{x)y' 



where I{x) stands for the integral 



I{x)= i'^^^du. 
uqo{u) 







Since the integrand clearly has no singularities at zero, the integral I{x) represents a smooth 
function of x. As it is easy to calculate, 

d{X, Y) _ ^ 
d{x,y) 

at the point x = y = 0, so that the transformation x ^ X, y Y is invertible on some 
neighborhood of the origin. Moreover, system (11), with the coefficients q2 and gs expressed 
in terms of go and qi by means of (12), in the new variables is reduced to the form 

fiX = -Y, nY = X, 

where the variable factor /x equals 1 at X = y = 0. Thus, the origin is a center of the system 
under consideration. 

As far as isochronicity is concerned, here it suffices to repeat what was said in this 
connection in the proof of the previous theorem: as it is easy to check, the field of the 
system commutes with the field (r, s) transversal to it in a neighborhood of the origin and 
defined by relations (|T^) . 

However, in this case isochronicity of the center follows also from other arguments. 
Namely, the above-mentioned factor n in the new variables X and Y has the form 1 — f{X)Y , 
where f{X) is the ratio qi{x)/X expressed in terms of X, and it is easy to verify by passage 
to polar coordinates that such a factor does not change the period of cycles of the harmonic 
oscillator. 

Closing the article, we return to the polynomial Abel systems and sum up some conclu- 
sions. First, as Theorem 7 and, especially. Theorem 8 demonstrate, such systems commute 
with transversal polynomial fields only in exceptional cases. Second, the very fact of com- 
mutation, as Theorem 9 or more general Theorem 10 asserts, guarantees the presence of 
a center at the origin and this center must be isochronous. Third, system of the form (^) 
and Theorems 3 and 4 devoted to them recall us that Theorem 8 describes far from all 
isochronous polynomial Abel systems. 

Problem 7. Find all polynomial Abel systems possessing isochronous centers. 
Acknowledgement We are thankfuU to S.Treskov and V.Vershinin for your help during 
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